The role of wire imperfections in micro magnetic traps for atoms 
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We present a quantitative study of roughness in the magnitude of the magnetic field produced by 
a current carrying microwire, i.e. in the trapping potential for paramagnetic atoms. We show that 
this potential roughness arises from deviations in the wire current flow due to geometric fluctuations 
of the edges of the wire : a measurement of the potential using cold trapped atoms agrees with 
the potential computed from the measurement of the wire edge roughness by a scanning electron 
microscope. 
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The use of micro or even nano-fabricated electrical de- 
vices to trap and manipulate cold atoms has attracted 
substantial interest, especially since the demonstration 
of Bose-Einstein condensate using such structures 0, Q • 
Compact and robust systems for producing EEC's, sin- 
gle mode waveguides and possibly atom interferometers 
can now be envisaged. The small size of the trapping ele- 
ments, usually current-carrying wires producing magnetic 
traps, and the proximity of the atoms to these elements 
(typically tens of microns), means that many complex 
and rapidly varying potentials can be designed Q. This 
approach has some disadvantages however. On the one 
hand it has been shown that atoms are sensitive to the 
magnetic fields generated by thermally fluctuatingcur- 
rents in a metal when they are very close 0, 0, IB IE Q • 
On the other hand, a time independent fragmentation of 
a cold atomic cloud has been observed when atoms are 
brought close to a current carrying micro- wire 0, |^ 0| . 

This fragmentation has been shown to be due to a po- 
tential roughness arising from distortions of the current 
flow in the wire 0,^3. It has also been demonstrated ex- 
perimentally that the effect of these distortions decreases 
with increasing distance from the wire |lOl [Tlj | . In an at- 
tempt to account for the observations, a theoretical sug- 
gestion has been made that the current distortions may 
be simply due to geometrical deformations, more specif- 
ically meanders, of the wire (IZj . In this paper we show 
that for at least one realization of a microfabricated mag- 
netic trap, using electroplating of gold, this suggestion is 
substantially correct. We have measured the longitudinal 
density variation of a fragmented thermal cloud of atoms 
trapped above a wire, and inferred the rough magnetic 
potential. We have also made scanning electron micro- 
scope images of the same wire and measured the profile of 
the wire edges over the region explored by the atoms. The 
magnetic potential as a function of position deduced from 
the edge measurements is in good quantitative agreement 
with that inferred from the atomic density. We suspect 
that this result is not unique to our sample or fabrication 
process and we emphasize the quantitative criterion for 
the necessary wire quality to be used for atom manipula- 
tion. 

The wires we use are produced using standard micro- 



electronic techniques. A silicon wafer is first covered by 
a 200 nm silicon dioxide layer. Next, layers of titanium 
(20 nm) and gold (200 nm) are evaporated. The wire pat- 
tern is imprinted on a 6 fim thick photoresist using optical 
lithography. Gold is electroplated between the resist walls 
using the first gold layer as an electrode. After remov- 
ing the photoresist and the first gold and titanium layers, 
we obtain electroplated wires of thickness uq — 4.5 /im 
with a rectangular transverse profile (see Fig. . A pla- 
narizing dielectric layer (BCB, a benzocyclobutene-based 
polymer) is deposited to cover the central region of the 
chip. On top of the BCB, a 200 nm gold layer is evapo- 
rated to be used as an optical mirror for light at 780 nm. 
The distance between the center of the wire and the mir- 
ror layer has been measured to be 14(1) /im. 

The magnetic trap is produced by a current / flowing 
through a Z-shaped microwire |0| together with an ex- 
ternal uniform magnetic field Bq (along the y-axis, see 
Fig. ^ parallel to the chip surface and perpendicular to 
the central part of the wire. The central part of the 
Z-wire is 50 ^m wide and 2800 ^m long. Cold *^Rb 
atoms, collected in a surface magneto-optical trap, are 
loaded into the magnetic trap after a stage of optical mo- 
lasses and optical pumping to the |F = 2, m = 2) hyper- 
fine state. The trap is then compressed so that efficient 
forced evaporative cooling can be applied. Finally, the 
trap is decompressed. Final values of / and Bq vary from 
200 mA to 300 mA and from 3 G to 14 G respectively 
so that the height of the magnetic trap above the wire 
ranges from 33 /xm to 170 /im. An external longitudinal 
magnetic field of a few Gauss aligned along z is added 
to limit the strength of the transverse confinement and 
to avoid spin flip losses induced by technical noise. For 
these parameters, the trap is highly elongated along the 
z-axis. The transverse oscillation frequency is typically 
oj±/{2tt) — 3.5 kHz and 120 Hz for traps at 33/j,m and 
170 /im from the wire respectively. 

The potential roughness is deduced from measurements 
of the longitudinal density distribution of cold trapped 
atoms. The atomic density is probed using absorption 
imaging after the atoms have been released from the final 
trap by switching off the current in the Z-wire (switching 
time smaller than 100 /is). The probe beam is refiected 
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FIG. 1: (a) : Z-wire used to produce the magnetic trap, (b) : 
Cross section of the wire in the (xy) plane. The wire is covered 
with a layer of BCB polymer, and a thin gold layer acting as 
a mirror. The origin of the coordinate system is taken at the 
center of the wire. 



by the chip at 45° allowing us to have two images of 
the cloud on the same picture. From images taken just 
after (500 fis) switching off the Z-wire current we infer 
the longitudinal density n{z) = JJ dxdyn{x , y , z) . We 
also deduce the height of the atoms above the mirror 
layer from the distance between the two images. The 
temperature of the atoms is determined by measuring the 
expansion of the cloud in the transverse direction after 
longer times of flight (1 to 5 ms). 

To infer the longitudinal potential experienced by the 
atoms, we assume the potential is given by 



V{x,y,z) = "14 (z) + Vhe,rm{x,y) 



(1) 



where T4iarm(2;, y) is a transverse harmonic potential. Un- 
der this separability assumption, the longitudinal poten- 
tial is directly obtained from the measured longitudi- 
nal density of a cloud at thermal equilibrium using the 
Boltzmann law Vz{z) = —k-QT\Ti[n{zy). To maximize 
the sensitivity to the longitudinal potential variations, 
we choose a temperature of the same order as the vari- 
ations ( r ~ 0.4 jiK for traps at 170 /im from the wire 
and T ~ 2.2 /iK for traps at 33 /xm from the wire). The 
separability assumption has been checked experimentally 
by deducing a z-dependent oscillation frequency from the 
rms width of the transverse atomic density at different 
positions. At 33 /im from the wire, there is no evidence 
of a varying oscillation frequency. At a height of 170 /^m 
above the chip, over a longitudinal extent of 450 /xm, we 
deduce a variation of the transverse oscillation frequency 
of about 13%. In this case, the assumption of separability 
introduces an error of Q.2k^T in the deduced potential. 

The potential experienced by the atoms is V = ^b\B\. 
To a very good approximation, the magnetic field at the 
minimum of Vharm is along the z-axis (for our parameters, 
the deviation from the z-axis is computed to be always 
smaller than 1 mrad) so that the longitudinal potential 
is given by Vz{z) — ^bBz{z). For a perfect Z-wire, the 
longitudinal potential is solely due to the arms of the wire 
and has a smooth shape. However, we observe a rough 
potential which is a signature for the presence of an addi- 
tional spatially fluctuating longitudinal magnetic field. A 



spatially fluctuating transverse magnetic field of similar 
amplitude would give rise to transverse displacement of 
the potential which is undetectable with our imaging res- 
olution and small enough to leave our analysis unchanged. 

In the following, we present the data analysis which en- 
ables us to extract the longitudinal potential roughness. 
In order to have a large statistical sample and to gain 
access to low spatial frequencies, one must measure the 
potential roughness over a large fraction of the central 
wire. In our experiment, however, the longitudinal con- 
finement produced by the arms of the Z-wire itself is too 
strong to enable the atomic cloud to spread over the full 
extent of the central wire. To circumvent this difficulty, 
we add an adjustable longitudinal gradient of which 
shifts the atomic cloud along the central wire. We then 
measure the potential above different zones of the central 
wire. We typically use four different spatial zones which 
overlap each other by about 200 /im. We then reconstruct 
the potential over the total explored region by subtract- 
ing gradients from the potentials obtained in the different 
zones. Those gradients are chosen in order to minimize 
the difference between the potentials in the regions where 
they overlap. 

We are interested in those potential variations which 
differ from the smooth confining potential due to the arms 
of the wire. We thus subtract the expected confining po- 
tential of an ideal wire from the reconstructed potential. 
To find the expected potential, we model the arms of the 
Z by two infinitesimally thin, semi-infinite wires of width 
250 ^m, separated by a distance I and assume a uniform 
current distribution. We fit each reconstructed potential 
to the sum of the result of the model and a gradient, us- 
ing the gradient, I and the distance h above the wire as 
fitting parameters. The fitted values of I differ by a few 
percent from the nominal value 2.8 mm, while we find 
h — 13(1) /im+d where d is the distance from the mirror 
as measured in the trap images. This result is consistent 
with the measured 14 /im thickness of the BCB layer. 

The potential which remains after the above subtrac- 
tion procedure is plotted for different heights in Fig. |21 
For a fixed trap height h (fixed ratio I/Bq), we have 
checked that the potential is proportional to the cur- 
rent in the wire; therefore we normalize all measure- 
ments to the wire current. The most obvious observa- 
tion is that the amplitude of the roughness decreases as 
one gets further away from the wire. The spectral den- 
sity of the potential roughness is shown in Fig. |21for two 
different heights above the wire. We observe that the 
spectrum gets narrower as the distance from the wire 
increases (see inset in Fig. (S)). This is expected since 
fluctuations of wavelength much smaller than the height 
above the wire are averaged to zero. At high wave vectors 
(fc > 0.07 ^m-i at 33 /im and k > OMfim^^ at 80 /im) 
the spectrum exhibits plateaus which we interpret as in- 
strumental noise. We expect this noise level to depend on 
our experimental parameters such as temperature, cur- 
rent and atomic density. Qualitatively, smaller atom- wire 
distances, which are analyzed with higher temperatures 
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FIG. 2; Rough potentials normalized to the current in the 
Z-shaped wire for different heights from the wire. Solid lines : 
potentials measured using cold atomic clouds. Dashed lines : 
potentials calculated from the measured geometric roughness 
of the edges of the wire. The different curves have been shifted 
by 6/iK/A from each other. 
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FIG. 3; Potential spectral density J„ — 

l/(27r/2) /^(V;(z)14(« + u))e^^''du at 33 ^tm from the 
wire (fat lines) and at 80 /xm from the wire (thin lines). Solid 
lines : potential measured using cold atomic clouds, the inset 
shows the curves on a linear scale. Dashed lines : potentials 
calculated from the measured geometric ffuctuations of the 
edges of the wire. These estimations of the spectral density 
are made with the Welch algorithm [T^ using windows half 
the size of the total explored region (1.6mm). 



should result in higher plateaus. This is consistent with 
the observation. 

In the following we evaluate the rough potential due to 
edge fluctuations of the central wire. For this purpose, 
the edges of the wire are imaged using a scanning elec- 
tron microscope (SEM), after removal of the BCB layer 
by reactive ion etching. Figure ^a) indicates that the 
function /, which gives the deviation of the position of 
the wire edge from the mean position y — ±Wo/2, is 
roughly independent of x. We make the approximation 
that / depends only on z. We deduce / from SEM im- 
ages taken from above the wire (see Fig.QJb)). To resolve 
/, whose rms amplitude is only 0.2 /im, we use fields of 
view as small as 50 /xm. The function / is reconstructed 
over the entire length of the central wire using many im- 
ages having about 18 /xm overlaps. As shown in Fig.0Jc), 
several length scales appear in the spectrum. There are 
small fluctuations of correlation length of about 100 nm 
and, more importantly, fluctuations of larger wavelength 
(60 to 1000 Aim). 

The geometric fluctuations of the edges of the wire in- 
duce a distortion of the current flow which produces a lon- 
gitudinal magnetic field roughness responsible for a po- 
tential roughness. To compute the current density in the 
wire, we assume a uniform resistivity inside the wire. We 
also assume that fluctuations are small enough to make 
the current density distortion linear in ]h|R^ where jhjR 
are the fluctuations of the left and right edge of the wire 
respectively. The current density is in the yz-plane and, 
because the rough potential is proportional to the longi- 
tudinal magnetic field, we are only interested in its y com- 
ponent, jy. Because of symmetry, only the part of jy(z, y) 
which is even in y contributes to in the a:2;-plane. 
Thus, only the symmetric component = l/2(fL + Jr) 



is considered. The Fourier component f^{k) of /+ in- 
duces a transverse current density |12| | 



j + {k,y)^ikf+{k) 



cosh(fcy) 



Wo Wo cosh(A:Wo/2) ' 



(2) 



As the distances from the wire we consider (33 to 176 ^va) 
are much larger than the thickness of the wire uq = 
4.5 /im we will in the following assume an infinitely flat 
wire. To efficiently compute the longitudinal magnetic 
field produced by these current distortions, we use the ex- 
pansion on the modified Bessel functions of second kind 
Kn{kx), which is valid for x > Wq/2. This expansion is, 
in the a;z-plane, 

B^{k,x) = - ^k{c2n{k) + C2n+2{k))K2n+l{kx) (3) 



where 



n>0 



C2„(fc) = {-ir—uo 

IT 



Wo/2 



hniky) j+{k,y) dy, (4) 



/„ being the modified Bessel function of the first kind. 
For small wave numbers k such that kWo -C 1, one ex- 
pects the Cn coefficients to decrease rapidly with n. In- 
deed, for kr <^ 1, In{kr) ~ (fcr)"/(2"n!). In our data 
analysis, only k wave vectors smaller than 0.07/im^^ are 
considered for which kWo/2 < 0.8 so we expect the c„ co- 
efficients to decrease rapidly with n. In the calculations, 
only the terms up to rt = 20 are used. 

Eauations l2l4l andl51then allow us to compute the fluc- 
tuating potential from the measured function /. In Fig.|21 
we plot the spectral density of the potential roughness 
calculated from / for two different heights above the wire 
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FIG. 4: Imperfections of the edges of the wire. In the SEM 
image of the wire taken from the side (a), one can see that 
the edge deviation function / is roughly independent of x. (c) 
Spectral density of / extracted from SEM images taken from 
the top as in Fig.(b). 

(33 and 80 /im) and compare it with those obtained from 
the potential measured with the atoms. For both heights, 
and for wave vectors small enough so that the measure- 
ments made with the atoms are not limited by experi- 
mental noise, the two curves are in good agreement. As 
we have measured the / function on the whole region 
explored by the atoms, we can compute directly the ex- 
pected potential roughness. In Fig.|21 this calculated po- 
tential roughness is compared with the roughness mea- 
sured with the atoms for different heights above the wire. 
Remarkably the potential computed from the wire edges 
and the one deduced from the atomic distributions have 
not only consistent spectra but present well correlated 
profiles. We thus conclude that the potential roughness 
is due to the geometric fluctuations of the edges of the 



wire. The good agreement between the curves also vali- 
dates the assumption of uniform conductivity inside the 
wire used to compute the current distortion flow. 

In conclusion, we have shown that the potential rough- 
ness we observe can be attributed to the geometric fluc- 
tuations of the wire edges. Fluctuations at low wave vec- 
tors, responsible for most of the potential roughness, cor- 
respond to wire edge fluctuations of very small amplitude 
compared to their correlation length. We emphasize that 
a quantitative evaluation of these wire roughness compo- 
nents demands dedicated measurement methods. 

Furthermore, wire edge fluctuations put a lower limit 
on the possibility of down-scaling atom chips. For a given 
fabrication technology the wire edge fluctuations are ex- 
pected to be independent of the wire width Wq. Thus as- 
suming a white noise spectrum, the normalized potential 

roughness Kms/^ varies as l/W^' for fixed ratio /i/Wq, 
h being the distance to the wire In order to reduce 
the potential roughness, one must pay careful attention 
to edge fluctuations when choosing a fabrication process. 
For example, we are currently investigating electron beam 
lithography followed by gold evaporation. Preliminary 
measurements indicate a reduction of the spectral den- 
sity of the wire edge fluctuations by at least two orders 
of magnitude for wave vectors ranging from 0.1 /im"^ to 
10 /im~^ . This should allow us to reduce the spectral den- 
sity of the potential roughness by the same factor unless a 
new as yet unobserved phenomenon such as bulk inhomo- 
geneity sets a new limit on atom chip down-scaling |l5j |. 
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